Study of B c J/ipK decays in the perturbative QCD approach 
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In this note, we calculate the branching ratio of B c — > J/ipK in the framework of perturbative 
QCD approach based on kr factorization. This decay can occur only via tree level diagrams in the 
Standard Model. We find that the branching ratio of B c — > J/ipK is about (1 — 3) x 10 -4 . The 
large branching ratio and the clear signals of the final states make the measurement of B c — > J/ipK 
easily at LHC-b experiments. 
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The B meson rare decays provide a good place for testing the Standard Model (SM), studying CP violation and 
looking for possible new physics beyond the SM. The theoretical studies of B u ,d mesons decays have been studied 
widely in the literature, which are strongly supported by the detectors at the e + e~ colliders, such as the CLEO, 
BaBar and Belle. At the era of the Large Hadron Collider(LHC), there is still a room for B physics. LHC beauty 
experiments (LHCb) will extend the B-physics results from the B factories by investigating decays of heavier B hadrons, 
such as B s and B c mesons. It is estimated that about 5 x 10 10 B c mesons can be produced per year at LHCp], 0], 
so the studies of B c meson rare decays are necessary in the next a few years, it will highlight the advantages of B 
physics. 

In this paper, we study the rare decays B c — > J/ifjK in the Perturbative QCD approach (PQCD) [3]. In SM, 
B c — > J/ipK decays occur through only the tree level diagrams and so there is no CP violation in this channel. 
Therefore, the measurements of this channel might provide a ground for investigating new physics effects. 
I 1 For the decay B c — > J/ipK, the related effective Hamiltonian is given by KM 
Oh' 

H cS = -^V us V* b {Ci(p)OiGi) + C 2 (/i)0 2 ( M )} + H.c, (1) 
where Ci (/x) are Wilson coefficients at the renormalization scale fi and Oi are the local four-quark operators 

> ; 

Ox = (biCj)v-A(UjSi) V -A, 2 = (biCi) V -A{u Sj)v-A- (2) 
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Amplitude- Jd 4 k 1( l 4 k 2 d 4 k 3 Tr[C(t)$B c (fci)^(fc 2 )$K(fc 3 )ff(fci,fc 2 ,fc3,0]e" S(t) . (3) 

> : 

In our following calculations, the Wilson coefficient C(t), Sudakov factor Si(t)(i — B c ,J/ip,K) and the non- 
perturbative but universal wave function can be found in the Refs. [BI-[lol|. The hard part H are channel dependent 
i but fortunately perturbative calculable, which will be shown below. 

Fig. [T] shows the lowest order Feynman diagrams to be calculated in the PQCD approach where (a) and (b) are 
factorizable topology, (c) and (d) are nonfactorizable topology. After a straightforward calculation using the PQCD 
formalism Eq.([3]), we obtain the sum contributions of (a) and (b) 

M a [C] = -^ttCVM1 c f dx 2 f b 1 db 1 b 2 db 2 x{[r b r(l-r 2 )^,(x 2 ,b 2 )~2r(l-r 2 )x 2 ^(x 2 ,b 2 ) 



Here i and j are SU (3) color indices. Then the calculation of decay amplitude is to evaluate the hadronic matrix 
elements of the local operators. 

In the PQCD approach, the decay amplitude can be written as: 



>o Jo 



+ (l-r 2 )x 2 ^(x 2l b 2 ) - 2r b {l-r 2 )^(x 2 M)]as(tl)h£\x 2 MM)cM-S B Atl) - Sj/^JP^) 
+ [r 2 (l - r 2 )(r c - l)^{x 2 ,b 2 ) ~ r c (l - r 2 )^(x 2 , b 2 )]a s {tl)hf\x 2 , h, b 2 ) 

xeM-SBXtl)-S.j,M)\C(ti)}, (4) 
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(d) 



FIG. 1: The lowest order diagrams for B c — > J/ipK decay. 



where r = M Jf ^/M Bc ,r b = M b /M Bc ,r c = M c /M Bc and 

t\ = max(M BV /(l - x 2 )(r c - r 2 ), M B ^\r 2 - (1 - r 2 )x 2 |, l/b u l/b 2 ), 



The functions 



t\ = max(M B y/r* +r c - r 2 , 1/fei, l/6 2 ). (5) 

/ii 1) (a; 2 ,6 1 ,6 2 ) = 5 t (ar 2 )if (M BA /(l - ^)(r c - r 2 )M 
x [0(6a - - (1 - r 2 )x 2 )I (M B ^r 2 b - (1 - r^h) 

x K (M BV /r2-(l-r 2 )a;2b2) + (&i <-> 62)], (6) 
/ii 2) (^2,6i,6 2 ) = S t (r c )K (M B y/(l - x 2 )(r c - r 2 )b 2 ) 



x [9(b 2 - b 1 )I Q (M B ^r 2 +r c - ?5i)Jf (M BV /^ + r c - r 2 fe 2 ) + (61 6 2 )] , (7) 

come from the Fourier transformation of propagators of virtual quark and gluon in the hard part calculations. For 
the non-factorizable diagrams (c) and (d), all three meson wave functions are involved. Their total contribution is: 

16 f 1 f°° 

M C [C] = -— ir C F fB e M% e / dx 2 dx 3 / b 2 db 2 b 3 db 3 cj ) £(x 3 ) 

c J " 

x{[(l - r 2 )(l - r 2 - r c - x 3 + r 2 r c + r 2 x 3 )^(x 2 ,b 2 ) + r(l - r 2 )(l - r c - x 2 )(t>%{x 2 ,b 2 )\ 
a s (^)/ l ( 1 )(x 2 ,x 3 ,6 2 ,6 3 )exph5 Bc (^) - Sj,^{t\) - S K {t^]C{t\) + [(r 2 - 1)(1 + r 2 - 2r c - x 2 + x 3 
-r 2 x 2 - r 2 x 3 )cj)^{x 2 ,b 2 ) + r(l ~ 7 " 2 )( 1 ~ r c - x 2 )$j ) (x 2 , 6 2 )]a s (t 2 )/i^ 2) (x 2 , x 3 , & 2 , & 3 ) 
exp[-S Bc (t 2 ) - 5. w (t 2 ) - S^(t 2 )]C(t 2 )}, (8) 

where 



= max ( M B ^J\F^\, M B ^{r c - r 2 )(l - x 2 ), l/b 2 , 1/63), 
t 2 = max(M B A /|F, 2 2 J, MbV^ - r 2 )(l - x 2 ), I/62, 1/63). ( 9 ) 



^ ) (a;2,a; 3 ,6 2 ,6 3 ) = 



0(b 2 - b 3 )I (M B ^(l - x 2 )(r c - r 2 )b 3 )K Q (M B ^(l - x 2 )(r c - r 2 )b 2 ) 

} ( K (M B F (j) b 3 ), for F 2 } > 0\ 
+ {b2 ~ b3) f X {^HM B j\F^b 3 ), for F 2 )<0 J' ^ 
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and Fq) 's are denned by 



^( 2 i) = (1 - r c - x 3 + r 2 x 3 )(x 2 - 1), 
Ff 2) = (1 - x 2 )(r c - r 2 - x 3 + r 2 x 3 ). 



The total decay amplitude is then 
and the decay width is expressed as 



CiM c + (-d + C 2 )M a f K 



(11) 



(12) 



j/^ + )| 



(13) 



The following parameters have been used in our numerical calculation 

M Bc = 6.286 ± 0.005GeV, M b = 5.2GeV, M J/V; = 3.097GeV,M c = 1.82GeV, 
/ Bc = 489 ± 4MeV, f K = 160MeV, f J/4 , = 405 ± 14MeV, 
t Bc = (0.46 ±0.07) x 10~ 12 s, |K S | =0.2257 ±0.0021,1^1 = (41.6 ±0.6) x 10~ 3 . 



(14) 



If not specified, we shall take their central values as the default input. We have taken the constituent quark masses 
M b = 5.2GeV and M c = 1.82GeV from the ISGW2 modelfT^. As noted in Ref.[H|], this choice for M b and M c 
satisfies approximately the relation Mf, = M c + 3.4GeV, which is consistent with the well known formula relating 
the pole masses M^^oie and M CiPO / e in the Heavy Quark Effective Theory. We also update some parameters used in 
previous works by taking the values from the latest Particle Data Group publication and lattice QCD simulations. 

Our numerical analysis shows that \C\M c j ' {{\Ci+C2)M a fk)\ = 17%, which means that the dominated contributions 
to the branching ratio of B c — > J/tpK decays come from the factorizable topology [(a) and (b) in FigJT]. We list our 
numerical results on the branching ratio for B c — > J/ipK decays in Table HI where ui is the parameter in the wave 
function of the J/tp meson(see Ref.Q for details). From the numbers in Table U one can find that the branching ratio 
of B c — > J/ipK decays is sensitive to the parameter u>. The branching ratio of B c — > J/tpK decays is more sensitive 
to the quark masses, especially the c-quark's mass, as shown in Tabic [Til Therefore the B c — > J/ipK decays provide 
a good platform to understand the wave function of the J/ip meson and the constituent quark model. Besides the 
uncertainty from the parameter ui and the quark masses, we find the uncertainty of the decay constant /;u will bring 
about 7% uncertainty to the branching ratio of B c — > J/tpK decays. We investigate the branching ratio's dependence 
on the hard scale t in Eq. ([3]), which characterize the size of next-to- leading order contribution. The branching ratio 
is shown in Table [TTT1 with those uncertainties. By changing the hard scale t from 0.8t to 1.2£, we find the branching 
ratio for B c — > J/ i/jK decays changes little as shown in Table IIII1 that mean the uncertainties in the next-to-leading 
order contributions can be neglected for this decay mode. The value of Aqcd also affect the branching ratio of 
B c —> J/ipK decays, we have taken Aqcd = 250MeV at Nt = 4 as our default input. A recent determination of 
Aqcd gives A™^* = 234 ± 26MeV[17\. Our result is Br(B c -> J/ipK) = 1.70 x 10~ 4 , if Aqcd = 208MeV; and 
Br(B c — > J/rpK) = 2.02 x 10 -4 , if Aqcd — 260MeV. Considering the uncertainties from the input parameters, our 
results on the branching ratio of B c —> J/ipK decays are generally in the range of (1 — 3) x 10~ 4 . The large branching 
ratio and the clear signals of the final states make the measurement of B c — > J/i/jK easily at the LHC-b experiments. 





uj = 0.5GeV 


id 


= 0.6GeV 


0J 


= OJGeV 


Br (Be - 


-¥ J/ll)K + ) 


1.52 


1.96 


2.44 



TABLE I: Branch ratios in the unit 10 4 for different uj 

In the literature, there already exist a lot of studies on B c — > J/ifjK decays [l6l Il8l - [2"4l |. we show their results in 
Table HVl All the previous works on this decay model were based on naive factorization, which, as exp ected, can be 
quite accurate for the B c meson, since the quark-gluon sea is suppressed in the heavy quarkonium[20j. In the naive 
factorization approach, the decay amplitude can be expressed in terms of the hadronic transition form factors and of 
the leptonic decay constants. For the B c — > J/tp transition form factors, there is a large difference among the previous 
worksQ. The authors of Ref.pjl [2(| calculated the form factors in the frame work of QCD sum rules, as argued by 
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M c = 1.72GeV 


M c = 1.82GeV 


M c 


= 1.92GeV 


Br(B+ - 


-> J/ipK+) 


3.80 


1.96 


1.19 



TABLE II: Branch ratios in the unit 10" 4 for different M c using uj = 0.6GeV and M b = 5.2GeV 



Scale 


Br(B+ -> J/ijjK + ) 


0.8 t 


2.10 x 10~ 4 


t 


1.96 x 10~ 4 


1.2 t 


1.91 x 10~ 4 



TABLE III: Branching ratio using cj = O.GGeV with uncertainties 



the authors, the sum rule estimates of form factors are taken at zero transfer squared, while the dependence on q 2 
is beyond the reliable accuracy of the method. The authors of Ref . (15, l22l-[24| presented their calculations based on 
the constituent quark model, where, the B c decay form factors can be expressed through the overlap integrals of the 
meson wave functions. In Ref.[2l[, Chang and Chen used a very complicated approach to evaluate the weak current 
matrix elements. It is expected that the approach is available as long as the mesons in the initial and final states 
are of weak binding, for example, for the decay B c — > J/ipK. After updating the model parameters, Ivanov et al. 
studied exclusive nonleptonic and semileptonic decays of the B c meson within a relativistic constituent quark model 
in Ref. fl^ . For the branching ratio of B c — » J/ipK decays, their results are close to our predictions. In general, the 
results of the various model calculations on B c — > J/ip decays are of the same order of magnitude. Eventually, the 
variations among the theoretical predictions should be discerned by experimental results of the B c decays, where, one 
can learn much about the decay mechanism. 
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